Manifestly covariant classical correlation dynamics II. 
Transport equations and Hakim equilibrium conjecture 



OS . 

o ■ 
o ' 

<n : 

>v 



os 



i 

c 

o 
o 



> 

On 

in 
o 

OS 

o 



x 



Chushun Tian 

Institut fur Theoretische Physik, Ziilpicher Str. 77, Koln, D-50937, Germany 

(Dated: May 29, 2009) 

This is the second of a series of papers on the special relativistic classical statistical mechan- 
ics. Employing the general theory developed in the first paper we rigorously derive the relativistic 
Vlasov, Landau and Boltzmann equation. The latter two advocate the Jiittner distribution as the 
equilibrium distribution. We thus, at the full microscopic level, provide a support for the recent 
numerical finding [D. Cubero et. at, Phys. Rev. Lett. 99, 170601 (2007)] of the special relativis- 
tic generalization of the Maxwell-Boltzmann distribution. Furthermore, the present theory allows 
us to rigorously calculate various correlation functions at the relativistic many-body equilibrium. 
Therefore, the relativistic many-body equilibrium conjecture of Hakim is justified. 

PACS numbers: 03.30.+p, 52.25.Dg 



I. INTRODUCTION 

The kinetic theory is one of the pillars of studies 
of relativistic transport phenomena in various systems 
ranging from star clusters or galaxies [H, 0, 0, 13 to 
plasmas in fusion quantum chromodynamics [a, 0] 
and graphene discovered very recently [8|. The mani- 
festly covariant counterparts of various classical trans- 
port equations, such as the kinetic equation of Vlasov, 
Fokker-Planck, Landau and Boltzmann, were proposed 
long time ago (for a review see, for example, Ref. [§]), 
and their applications nowadays have been well docu- 
mented [Tol . These manifestly covariant transport 
equations have received justifications from various mi- 
croscopic approaches 0, H, H, HI, [H, EH EH], ft is 
rather typical in rederiving relativistic transport equa- 
tions (for examples, see [H EU, EU) that the Liouvil- 
lian dynamics of complete many-body distribution func- 
tion is bypassed and, crucially, the truncation approx- 
imation is resorted to. Therefore, despite of the great 
success of this kind of microscopic approaches 0, [Tj, EH] 
a fundamentally important problem remains unsolved. 
That is, will a manifestly covariant transport equation be 
compatible with the Liouvillian dynamics of the complete 
many-body distribution function? This issue has been ad- 
dressed continuously for several decades by many workers 
@, i, i, [H GJ, E3, El, d and been far less advanced. 
This is by no means of pure theoretical interests, rather, 
may find considerable practical applications. Indeed, ex- 
periences in the Newtonian physics have shown that to go 
beyond the weak coupling and Markovian approximation 
is inevitable in order to understand collective dielectric 
effects in electromagnetic plasmas [2(| and the infrared 
divergence of the Fokker-Planck equation describing the 
stellar dynamics [2ll |22| . In an insightful paper Kandrup 
first realized that a Liouville equation of the complete 
many-body distribution function admits an exact closed 
(nonlinear) kinetic equation which is satisfied by the re- 
duced one-body distribution function and manifestly co- 
variant Q. Various special transport equations, remark- 
ably, may be unified within this general kinetic equation. 



Unfortunately, he encountered principal difficulties when 
proceeded further to derive the general kinetic equation 
explicitly. Such a big gap is filled in the first of this 
series of papers (denoted as Paper I) by the developed 
manifestly covariant classical correlation dynamics [23j |. 
One of the main subjects of the present paper is to re- 
cover various special transport equations systematically 
from the exact, manifestly covariant and closed kinetic 
equation [24j] found in Paper I, which we call the general 
transport equation. 

The existence of the general transport equation has im- 
portant implications on several aspects of classical special 
relativistic (nonequilibrium) statistical mechanics which 
are undergoing intensified studies. First of all, the math- 
ematical structure possessed by the general transport 
equation is in excellent agreement with various approxi- 
mate transport equations obtained by other microscopic 
approaches @, S El El Hi, HI, El S , but remarkably 
different from that proposed by Horwitz and coworkers 
|l7l ]. Recently triggered by the latter there have been 
many scientific activities [2a , |27| ] searching for the spe- 
cial relativistic generalization of the Maxwell-Boltzmann 
distribution. The general transport equation reinforces 
the concept well established by approximate transport 
equations. That is, the Jiittner distribution serves as 
the relativistic (one-body) equilibrium in dilute systems. 
Therefore, it is suggested that other alternatives to the 
Jiittner equilibrium [17J might be specific to the (deter- 
ministic) relativistic many-body dynamics. 

Secondly, in past years the relativistic Brownian mo- 
tion and diffusion have experienced considerable concep- 
tual developments and found important practical appli- 
cations [281 ] . Long time ago it was known that in the 
Minkowski spacetime nontrivial Lorentz invariant Marko- 
vian processes do not exist [H [3(| • It turns out that the 
relativistic Brownian motion is interpreted as relativistic 
Markovian processes in the fi phase space, and the latter 
is completely described by a relativistic Fokker- Planck- 
type equation. So far these observations have been in- 
vestigated thoroughly at the level of one-particle physics 
[III , and Paper I is the first to substantiate these impor- 
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tant observations at the level of the genuine relativistic 
many-body physics. Indeed, the general transport equa- 
tion arises from (i) the thermodynamic limit namely the 
particle number N — > +00 and, (ii) that at given (global) 
proper time particles lose the memory of the history of 
the entire system. The condition (ii) is in sharp con- 
trast to the deterministic relativistic many-body dynam- 
ics, where the particle interaction is profoundly nonlocal 
in spacetime, and serves as the many-body dynamical 
origin of the Markovian processes in the /j, phase space. 
In particular, in Paper I the proper time parametrized 
equation in the 8-dimensional /1 phase space is rigorously 
justified, which was first obtained by Hakim [30| . and im- 
portant roles of which have very recently been reinforced 
[3l|. From these perspectives the theory presented in 
this series of papers may be considered as a microscopic 
approach to (classical) special relativistic Brownian mo- 
tion complementary to the one based on the relativistic 
Langevin equation. In particular, it might be proven to 
be a useful technique in exploring the concepts such as 
relativistic noises, friction and fluctuation theorem. 

There is an adjacent important yet unsolved problem 
that may be explored in the present theoretical scope. 
That is, to formulate the relativistic many-body equi- 
librium. In the notable critical analysis [9( Hakim con- 
jectured that at equilibrium there might exist an infi- 
nite Lorentz invariant hierarchy of correlation functions 
which is invariant under the spacetime translation and is 
merely determined by the Jiittner distribution. Unfortu- 
nately, further progresses have been impeded by the trun- 
cation approximation intrinsic to various microscopic ap- 
proaches (for examples, Refs. 0, HH), and the (dis)proof 
so far has been missing. The general principles given 
in Paper I paves the way towards justifying this conjec- 
ture. There the hierarchy of (physical) correlation func- 
tions is found explicitly, which is merely determined by 
the (physical) one-body distribution function. Provided 
that the general transport equation admits an equilib- 
rium distribution, then the hierarchy of equilibrium cor- 
relation functions is uniquely determined. To carry out 
this program for rarificd electromagnetic plasmas consti- 
tutes another main subject of the present paper. 

The present paper is written in the self-contained man- 
ner. The readers that would not like to access to the 
mathematical foundation may skip Paper I. The paper 
is organized as follows: Sec. |TT] is an exposition of the 
main results of Paper I, and the exact starting point of 
the present paper is pointed out. The rest is devoted to 
applications in relativistic plasmas with electromagnetic 
interactions. In Sec. lIIII we derive the relativistic Vlasov, 
Landau and Boltzmann equation from the general trans- 
port equation found in Paper I. In Sec. II VI we show that 
the present theory fully agrees with the recent numer- 
ical finding [26| and advocates the Jiittner distribution 
as the special relativistic generalization of the Maxwell- 
Boltzmann distribution. The Hakim equilibrium conjec- 
ture is justified, and the two-body equilibrium correlation 
function is exactly calculated. We conclude this series of 



papers in Sec. [Vj Some technical details are given in 
Appendix MCl 

Finally we list some of the notations and conventions. 
We choose the unit system with the speed of light c = 1 . 
To distinguish from the Minkowski 4-vector we use the 
bold font to denote the vector in the Euclidean space. 
Greek indices running from to 3 are further used to 
denote the component of the former. The Einstein sum- 
mation convention is applied to these indices. The 4- 
dimensional Minkowski space is endowed with the met- 
ric rf v = diag(l, —1, —1, —1) . The scalar product of two 
4- vectors is defined as a ■ b = if* 'a^&„ = . In partic- 
ular, a ■ a = a 2 . In addition to the usual mathematical 
symbols we use the following notations: 

<9 M , covariant derivative: <9 M = d/dx^ 1 ; 

d 4 2 , volume element in 4-dimensional Minkowski 
space: d A z = dz°dz 1 dz 2 dz 3 ; 

d 3 z , volume clement in 3-dimensional Euclidean 
space: d 3 z = dz 1 dz 2 dz 3 ; 

S^ d \f) , e?-dimensional Dirac function: 

6(x) , Heaviside function ; 

d£ M , differential form of spacclikc 3-surfacc: 

dEp = ^e^upx dx v A dx p A dx x with e^px 
being ±1 when (nvp\) is an even (odd) 
permutation of (0123) and being otherwise ; 

Xi [<;] , world line of particle i ; 

Xj(^) , 4-position of particle i at proper time s . 



II. MAIN RESULTS OF GENERAL THEORY 

This section is devoted to presenting the exact start- 
ing point that underlies the entire analysis of the follow- 
ing sections. For this purpose we first briefly review the 
manifestly covariant correlation dynamics developed in 
Paper I, and introduce all the mathematical objects used 
throughout this paper. It should be stressed that at each 
step of the manipulations below the manifest covariance 
is preserved. 



Correlation dynamics of r-parametrized 
evolution 



Consider a microscopic system composed of N iden- 
tical (classical) particles with mass m , the dynamics of 
which is formulated within the action-at-a-distance for- 
malism [32j. More precisely, the history of the system 
is described by a bundle of N particle world lines which 
solve the following relativistic motion equations: 



dxl 

dn 
dn 



m ~ 1 ' 

N 



(1) 
(2) 



Here xf (t^) , uf(ri) ,pf(Tj) are the 4-position, the 4- 
velocity and the 4-momentum vector of particle i de- 
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pending on the proper time Tj , respectively, and F£ is 
the force acting on particle i by particle j . Notice that 
F^(xi,pi) functionally depends on the world line of par- 
ticle j namely Xj [Tj] . For simplicity we here consider the 
case where the external force is absent, and the interact- 
ing force F^. is conservative, i.e., 



d 

—Ft J (x l , Pl ) = Q. 

The mass-shell constraint: p\ — m 2 is preserved by 
Pi ■ Fij(xi,pi) = . 



(3) 



(4) 



Eqs. JT]) and @ suggest that in order to formulate 
a statistical theory of an ensemble of such systems the 
introduction of an 8-/V-dimensional T phase space and 
the associated probability (phase) density function T> is 
required. (The mass-shell constraint is absorbed into 
the distribution function.) Remarkably, this distribu- 
tion function differs from its Newtonian counterpart in 
that it is parametrized by N (rather than 1) proper 
times, and functionally depends on N particle world 
lines. For T> the probability conservation law gives 
N manifestly covariant Liouvillc equations. By fur- 
ther introducing an auxiliary "gauge" condition-to de- 
mand the N proper times to change uniformly-we obtain 
a manifestly covariant single-time Liouville equation of 
T>(xx,p 1 ,n + t, ■ ■ ■ ,x N ,p N ,T N + t;xi[^}, ■ ■ ■ jXjvfc]) : 



£ X> = , 



d_ 

&? 

where the LiouviUian £ is given by 
£ = £° + A£', 



(5) 



(6) 



1=1 



i<j 



A 4" = -{FtM'Pi^+Fjtixj'Pi)^ 



with £° and A£^- the free and the two-body interacting 
LiouviUian, respectively. Notice that here the dimension- 
less parameter A characterizes the interaction strength. 

For the evolution parametrized by r namely Eq. ([5]) 
the correlation dynamics analysis may be performed for 
a large class of realistic systems [23j . First, we define the 
following distribution vector: 



23 =({Z>i}, {%},•■• ,{Vn} = V).. 



(7) 



where the reduced s-body distribution function T> s is ob- 
tained by integrating out arbitrary (N — s) particle phase 
coordinates and because of this for each {T> s } there are 
Nl/[(N — s\)sl) components. With the help of this def- 
inition the BBGKY hierarchy is rewritten in a compact 
form: 



7^- £! V = ,! ' 



(8) 



Eqs. {TJ and ([8]) constitute the reduced distribution func- 
tion representation of the single-time Liouville equation 
([5|). To proceed further we introduce the so-called cor- 
relation pattern representation allowing a more delicate 
decomposition of the distribution vector or more gen- 
eral functions. A correlation pattern, denoted as |r s ) 
(or (r s |), describes the statistical correlation of given 
s-particle group. More precisely, consider an s-particle 
group (ii---i s ),s < N, then the correlation pattern 
is generally given by (P^l • • • |Pj) (or (P 1 |P 2 | • • • |Pj|), 
where Pi , • • • , Pj is a partition of (i\ • • ■ i 8 ) ■ It implies 
that in the statistical sense within Pj , 1 < i < j the par- 
ticles correlate with each other, while the particle groups 
Pi , 1 < i < j are independent. The reduced s-body 
distribution function in this presentation now reads out 
as T> s = J^r |r s ) (r s |2? s which is none but the cluster 
expansion. It should be stressed, however, that it dif- 
fers from a traditional one [3^, in that the distribu- 
tion functions depend on the particle world lines. This 
is, indeed, an important ingredient of the Klimontovich 
technique in the Newtonian context [l3| and was gener- 
alized to the special relativity-in a manifestly covariant 
manner-by Hakim Q. 

A particularly important correlation pattern is the so- 
called vacuum state: \T S ) = |0 S ) (or (r s | = (0 S |) , where 
the given s particles are (statistically) independent. In 
contrast, all the other correlation patterns are called cor- 
relation states. With this definition the vacuum and the 
correlation operator, denoted as V and C , respectively, 
are defined as follows: 



vir,) = £. 



o r r r |r r ) , 



C|r r ) = (1 - 5 0r r r ) \T r ) , (9) 



which project given functions onto the vacuum or the 
correlation state. 

Then, a scries of rigorous theorems establish the fol- 
lowing important properties. First of all, the distribution 
vector D is split into the kinetic component Ilk. 2) and 
the nonkinetic component fl n .k. 23 , i.e., 



23 =n k .23 +n n . k .23 



(10) 



and the latter is irrelevant for large global proper times. 
The former is further decomposed into the vacuum and 
the correlation state, i.e., 



n k .23 = vn k .23 + cn k .23 



(11) 



and the correlation state is fully determined by the vac- 
uum state. In the thermodynamic limit N — > +oo as 
written in the correlation pattern representation the evo- 
lution of VII k . 23 and CII k . 23 are expressed in terms 
of two infinite equation hierarchies. It is remarkable 
that both infinite equation hierarchies are determined 
merely by a reduced one-body distribution function 



4 



T>(x,p; Xr XtP \) which solves the exact closed equation: 



+ ufa M i - J d2Xl' 12 V(2;X 2 )^V(l;X 1 ) 
Y,j d2 - j ^'(l|V(r - £)V|1| • • • li) n T>{s; X s 



(12) 

Here we have introduced the shorthand notations: 
T>{i;Xi) = 'D(xi,pi\Xt Xi>Pi \) and di = d A Xid A pi , and 
the notation: Xt x ^ p \ stands for some world line passing 
through the phase point (x,p) . Notice that the two- 
body interacting Liouvillian XC' i j(xi,pi]Xj,pj) is a func- 
tional of the world lines X( Xip .j and X(^ x . p .\ . The oper- 
ator VTV determines the evolution of the vacuum state 
Vllk. 2) , and is given by the following functional equa- 
tion: 



VrV = V£V+ / daVfl(a)Vexp(-sVTV), (13) 
Jo 

f d7 a * 

V§(s)V= / — e- lzs V£(z)m°(z)C\Z'V , 
Jc 2lT 

where the contour C (in the complex plane) lies above 
all the singularities of the Laplace transform of T> , and 



£{z) = ^A™ +1 £'{C«n (z)£'}", (14) 



n=0 



i 



-iz - £° 



(15) 



To our best knowledge a kinetic equation similar to 
Eq. (|12p was first obtained by Hakim by using the weak 
coupling approximation [9]. A simplified equation of 
Hakim has very recently been found to play important 
roles in studies of the relativistic Brownian motion 13 ill. 



B. General transport equation and hierarchy of 
physical correlation functions 

The r-parametrized evolution may not be observable 
because T> is normalized in the 8iV-dimensional phase 
space. Rather, to match macroscopic observations the 
following physical distribution function: 



M{xi,pi, ■ ■ ■ ,x N ,p N ]Xi{<;} 



,x N 



N 

[<;)) = JUdnViw) 



is introduced which may be considered formally as 
the stationary solution to Eq. ([5]). As a re- 
sult, the nonkinctic component of J\f identically van- 
ishes. Moreover, in the thermodynamic limit A — > 
+oo both hierarchies (corresponding to the vacuum 
and the correlation state, respectively) is determined 
merely by a physical one-body distribution function 



f(x,p). Here by "physical" the normalization condi- 
tion: lini7v^+oo A -1 / S( g,(j4 dUfid^p f(x,p) = 1 is im- 
plied, where S is a spacelike 3-surface, and U 4 is the 
4-dimensional Minkowski momentum space. 

The stationary solution to Eq. (TT2"]) results in a gen- 
eral transport equation which is manifestly covariant and 
closed: 



£ 2 (g>[/! 



dH^p 2 u^\C' l2 f{2) /(l) = K[f] 



(17) 



where f(i) is the shorthand notation of f(xi,pi) , and the 
collision integral is given by 



K[/] = V" / dT, ll2 d i p2U l 2 l ■ ■ ■ j dT, llj d i pjU 



j>2 J V2®U$ 

x<l|V(T- £)V|1| ...|j> JJ/(i). 



(18) 



The solution to Eqs. (fT7)l and (JTSJ, in turn, uniquely de- 
termines the hierarchy of physical correlation functions. 
More precisely, given an arbitrary j-particle correlation 
pattern 1^(1, • • • ,j) ^ Oj , in the thermodynamic limit 
A — > +oo the physical correlation function (denoted as 

(rj lDToo) reads out as 

(iilC (19) 

„oo „ s oo n „ n 

= ds ds'J2 II / dH^piut TT/(fc) 

Jo Jo , JSigt/. 4 



n=j i=2+l ' 



k=l 



x<rj-|c£l (s - s')€(s')Vexp(-sVrV)|l| ■•• \n) , 

(The integration procedure: n£=j+i Js,®£/ 4 dY.^d^piu'f 
is defined as unity for n = j .) where 

^ (r) = I c £ e ~ tZTSk °^ < ^ = J c £ e ~ l2T ^) • 

(20) 

In practical applications a perturbative expansion with 
respect to A may be further performed for the colli- 
sion integral K[/] and the physical correlation function 
(1^ I D^oo . Such an expansion is represented by the di- 
agrams constructed out of the free propagator and the 
interaction vertex (see the inset of Fig. [1]). For the for- 
mer the matrix element reads out as 



(z'|e" T " a " |*> = M«' e- T < 9 »< 



(21) 



For the latter there are two types: For the first type (the 
middle in the inset of Fig. [I]) a particle joining the vertex 
from the right is annihilated (the dashed line), and the 
matrix element reads out as 



(z|AA Y |r 2 (i,j)) 

{Su'Sjf + 5ij>6ji>) 



dE^pjU^XC'i, (22) 
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FIG. 1: Diagrams representing the perturbative expansion 
of the collision integral: (a) the weak coupling approxima- 
tion and (b) the lowest order ring approximation. The circle 
stands for the physical one-body distribution function. In- 
set: free propagator (top) and two types of interaction vertex 
(middle and bottom). 



with \T2(i,j)) = or For the second type (the 
bottom in the inset of Fig. [J) no particles (joining the 
vertex from the right) are annihilated. "Switching on" 
the interaction (at the vertex) introduces the (statisti- 
cal) correlation, i.e., (£2(1, j)\ = irrespective of the 
"initial" correlation pattern (to the right of the vertex), 
i.e., |T2 («, j)} ■ The matrix element for this type of the 
interaction vertex reads out as 

(ij|A^; ,|r 2 (i,i)) = (6u>S jf + S^S^XC^ . (23) 

Eqs. (fTT)) - (flT))) justify the manifestly covariant Bogoli- 
ubov functional assumption [33l . l3~il . 13a ] . They constitute 
the complete set for describing various physical phenom- 
ena such as transport processes, (macroscopic) relativis- 
tic hydrodynamics, and (physical) correlations at equilib- 
rium, and serve as the exact starting point of subsequent 
sections. The remaining of this paper is, indeed, devoted 
to the applications of Eqs. (fT7| - (fll))) in classical relativis- 
tic plasmas with electromagnetic interactions. 



III. SPECIAL TRANSPORT EQUATIONS 

In this section we will consider a rarified electron 
plasma that is near the local equilibrium with the den- 
sity and the temperature in the local rest frame as po 
and T , respectively. We will justify that various (namely 



Vlasov 9], Landau [Hll] and Boltzmann [lEGIj]) man- 
ifestly covariant kinetic equations existing in literatures 
are unified within the general transport equation (|17[) . 
Particular attention will be paid to the additional ap- 
proximations made. In doing so we expect to clarify the 
context where they are applicable. 

It should be stressed that given a density po the present 
theory is applicable for moderate temperatures. For suffi- 
ciently high or low temperatures quantum statistics and 
QED processes dominate, and the complete treatment 
requires a quantum theory. Let us now estimate such a 
condition. On one hand, two classical electrons may ap- 
proach each other up to a distance ~ e 2 / (ksT) , where fcg 
is the Boltzmann constant. In order for a classical scat- 
tering theory to be applicable it must be larger than the 
thermal deBroglie wavelength h/y/mksT . Consequently, 
we find fc^T <C a 2 m with a = — |y the fine structure con- 
stant. On the other hand, the plasma become degenerate 
when the thermal deBroglie wavelength and the mean 

— 1/3 

distance ~ p are comparable. This implies that the 
quantum statistics can be ignored only for sufficiently 

—1/3 

high temperatures such that H/^mkBT -C Pq , i.e., 
ksT 3> h 2 p 2 / 3 /to . Thus, the temperature region for a 
classical theory to be applicable is 

^£°l <<kB T«a 2 m. (24) 

TO 

Notice that the inequality above imposes a restriction on 
the density, i.e., po <C (e 2 m/h 2 ) 3 . This implies that the 

fermi energy namely h 2 p^ 3 /m is the lowest energy scale 
and, in particular, is much smaller than the (classical) 

electromagnetic energy e 2 pQ . The inequality (j2~4"|) is 
the rigorous condition for the present classical theory to 
be applicable. 

A. General scheme 

Let us first prove an exact relation between the colli- 
sion integral K[/] and the physical two-body correlation 
function. For the latter setting (Tj\ to be (12| for the 
hierarchy (fll?]) we find 

<12|^oo (25) 

POO PS 00 71 P n 

= ds ds'Y] J] / dE^piuf T[f(k) 
Jo Jo n =ji=j+l J ^® U t k=l 

x(12|Cil (s - s')£(s')Vexp(-sVrV)|l| ■• • \n) . 
Notice that 

V(r-£)V= (26) 

pOO pS 

/ ds ds'VA£'c£l°(s-s')f(s') Vex P(- sVrv )- 
Jo Jo 

Inserting it into Eq. (fl8|) gives 

K[/]=/ dS M2 d 4 p 2 <(I|A£'|12)(12|^ co . (27) 
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This exact relation indicates that the partition of full dis- 
tribution functions in the correlation pattern representa- 
tion preserves the cluster expansion, which is guaranteed 
by the Clavin theorem in the Newtonian physics [37l ]. 

We remark that from Eq. (|27|) the collision integral is 
locally well defined provided that either interactions or 
(statistical) correlations are short-ranged in the space- 
like 3-surface passing through x\ and x 2 ■ A relativis- 
ts plasma with electromagnetic interactions, indeed, be- 
longs to the latter case. There, although the static 
transverse electromagnetic field is long-ranged, the De- 
bye screening of the longitudinal electromagnetic field, 
as we will show in Sec. IIV1 renders the correlation func- 
tion short-ranged with the correlation radius 



I k B T 
47re 2 po 



(28) 



which is the exact condition for the subsequent analysis 
to be applicable. 

Then, the low density limit introduces a substantial 
simplification of the collision integral K[/] . Indeed, be- 
cause the physical one-body distribution is proportional 
to the density in Eq. (jTTJ) the collision integral may be 
considered formally as a density expansion which, term 
by term, corresponds to the two-, three-body scattering 
and so on. In the low density limit all the higher order 
terms in this density expansion may be ignored. As a 
result, 



K[/] = / dZ,ad*p2U$ 



x(l|V(r-£)V|l|2)/(l)/(2). (30) 



As a result, in Eq. (|27p the integral over d£ M2 is domi- 
nated by the region around x\ of size Ad • Notice that 
the correlation radius Eq. ([28]) makes sense only if it is 
much larger than the mean distance between two nearest 

— 1/3 

electrons, which is order of p . This leads to a suf- 
ficiently small plasma parameter, i.e., e 2 pj/ 3 /ksT <C 1 • 
Combining with the inequality (124[) we obtain 



Furthermore, with the help of appropriate iteration for 
Eq. J2H|) we explicitly write down the operator VTV as 



v(r-£)v = £ vr [n] v. 



(31) 



n=l 



^2/3 



Po „ 2 1/3 i rp 2 

— — < e p ' < k B T < a m , 



(29) 



Here 



T[ n ] = / ds 2 / ds 4 --- ds 2n / ds 3 ds 5 --- ds 2n -i V5(s 2 )ViX(s 1 - s 2 - s 3 )V 

Jo Jo Jo Jo Jo Jo 

X0(s 4 )Vil(s 3 -84,- S 5 )V • • • fl(S2„)Vil(s 2 n-l ~ S 2n ~ S 2 „+l)V , (32) 

I 



where s% = s 2n _|_i = , and the operator Vil(s)V = 
Vexp{sV£V}. Notice that in order for the matrix 
element (r r |V§(s)V|r^.,) not to vanish the condition: 
r r = r ,?* > 2 and IV = r /,r' > 2 must be met. 
For IVi with n > 2 there are more than one particles an- 
nihilated, which leads to higher order density corrections 
as the interaction strength A is compensated by a den- 
sity factor associated with the annihilated particle. For 
this reason in the expansion of Eq. (f3"Tj) only the leading 
term (n = 1) is kept. Furthermore, because in the A- 
expansion of V1XV the higher order terms are associated 
with the particle annihilation resulting in higher-order 
density corrections, the replacement: iX(— s) — * it°(— s) 
may be made. Consequently, we obtain 

K[f] = f cE; l2 d 4 p 2 7^ (33) 

J^2®U^ 
r- OO 

x / d s (l|Vfl( s )v£l (~ s )V|l|2)/(l)/(2) 
Jo 



as the (formal) leading order density expansion of the 
collision integral (|3"D|) concerned. 

The simplified collision integral l|33p can be expressed 
in terms of the A-expansion. It is important that for 
this expansion the interaction strength A is not compen- 
sated by the density factor because no particles are anni- 
hilated. In the remaining part we show that keeping such 
an expansion up to the A/ A 2 term results in the mani- 
festly covariant Vlasov/Landau equation, while keeping 
the entire expansion results in the manifestly covariant 
Boltzmann equation. 

Finally let us present a summary of the approximations 
to be used in the subsequent analysis that implement the 
scheme outlined above. One is the so-called relativistic 
impulse approximation There, the phase trajectory 

X( XlP ) in Eqs. (fT7|) and (fTg)) is given by 

X( x ,p) = x[s] — x 11 - u^s , (34) 

where we choose the proper time origin to be the mo- 
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ment at which the world line passes through x with the 
given 4-velocity u M . The other is the traditional hydro- 
dynamic approximation assuming that the physical dis- 
tribution function f(x,p) varies over a spatial (tempo- 
ral) scale much larger than Ad {^p 1 = y/m/ '47re 2 po) ■ 
Eq. (|33[) . in combination with these two approximations, 
is the starting point of subsequent analysis. 

B. Mean field approximation: Relativistic Vlasov 
equation 

Let us start from the simplest case namely to keep the 
A-expansion up to the first order. That is, we neglect the 



collision integral. Consequently, we obtain 

| - A / d^ 2 d 4 P2 u%t' l2 f(2) 1 /(l) = 0(35) 

[ Jn 2 ®uj J 

According to the second term of Eq. (j35[) . the physical 
one-body distribution function is driven by the mean field 
formed by all the other particles. For this reason, to 
keep the interaction expansion up to the leading order is 
called mean field approximation. To fully determine the 
mean field we use the relativistic impulse approximation 
Eq. ((35). A s a result, 



■o 2 2 r jAu 

m J (27r) 4 



Q l2 (k) ee ± ^S(k-p 2 )[k^p 1 -p 2 )^p^k-p 1 )]^-S(k-p 1 )[k^ Pl -p 2 )-p^k-p 2 )]^Y (36) 
The derivation is given in Appendix IA 21 

C. Weak coupling approximation: Relativistic Landau equation 



In this part we will consider the so-called weak coupling approximation to Eq. (|33|) . That is, the A-expansion is kept 
up to the second order [Fig. Q] (a)]. The collision integral thereby obtained is denoted as ~K\[f] ■ 

1. Collision integral 
Under the weak coupling approximation Eq. ()33[) is simplified as 



DO 



Kj[/] = A 2 / da I d^ 2 d 4 p 2 u 2 'C[ 2 VUs)£'i2V 12 (-s)f(x2,P2)f(x 1 ,p 1 ), (37) 

JO J£ 2 ®U* 

where the propagator U^-(s) is defined as 

TJO- (s) = exp[- S «^ + up M )} . (38) 

As discussed in Sec. IIII Al because of the short-ranged two-body correlation the integration over dS M 2 is dominated 
by a region of size Ad ■ Applying the hydrodynamic approximation gives 



dY>^ 2 d p 2 u 2 L (• • ■)f(x 1 ,p 1 )f(x 2 ,p 2 ) w / (K^d p 2 u^ {• ■ ■)f{x u p 1 )f{xi,p 2 ) (39) 
for the collision integral Ki[/] . That is, the collision integral is local in spacetime: 

poo P 

K 1 [/] = A 2 / da di:^ 2 d i p 2 u^C\ 2 t\ 2 (a)C' 12 t\ 2 (-a)f(x u p 2 )f(x ll p 1 ). (40) 

Jo Jt. 2 ®u% 

With the relativistic impulse approximation further used we obtain 



dp: 



f{xi,p 2 )f(xuPi), (41) 



2e\u 1 -u 2 ) 2 I d 4 kS(k- Ul )5(k-u 2 )-^— ^ (42) 
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after straightforward but tedious calculations, which are detailed in Appendix [B] Eq. (|4Tj) is the relativistic Landau 
collision integral [3611. It was first justified by Klimontovich at the full microscopic level albeit in a nonmanifestly 
covariant manner [l3l [Hj]. Notice that although the relativistic Landau collision integral is divergent, it formally 
admits the Jiittner distribution as the unique (local) equilibrium distribution (see Sec. IIV Al for detailed analysis), 
and gives the relaxation time (up to a numerical factor) as y m^ksT) 3 / (e po) S> ^p 1 ■ 



2. Logarithmic divergence of the collision integral 

However, the collision integral Ki [/] suffers from both the infrared and the ultraviolet divergence. Indeed, with the 
integral over wave vector carried out Eq. (|42|) gives (see Appendix [C] for details) 

= -27re 4 | ^ ( Ul • u 2 f[{ Ul ■ u 2 f - \]^ 2 

x {[(ui ■ u 2 f ~ l}g^ + «< + «) - {ui • u 2 )«U2 + «)} • (43) 
I 



A similar logarithmic divergence of the collision integral 
was first noticed by Landau (3§] in the context of non- 
relativistic plasmas (the so-called Coulomb logarithm). 
Following the prescription of Landau, from the practical 
viewpoint, in order to describe nonequilibrium processes 
near (local) Jiittner equilibrium it suffices to substitute 
appropriate ultraviolet (infrared) cutoff /c max (fcmin) into 
the collision integral since the divergence is logarithmic. 
To further estimate these cutoffs in the relativistic con- 
text we notice that is the minimal distance as two 
classical electrons approach each other. At such a dis- 
tance the kinetic energy and the interaction becomes 
comparable, i.e., 

fc max ~ ■ (44) 

For moderate temperatures fcsT <C a 2 m to heal the 
ultraviolet divergence is well within the reach of the 
present theory. Physically, the weak coupling approxi- 
mation namely the leading order A-expansion of Eq. (|33|) 
accounts for the small angle scattering, but fails in de- 
scribing large angle scattering. To implement this one 
needs to sum up the entire A-expansion of Eq. Q33p or the 
diagrams shown in Fig. O (Let us keep in mind that the 
collision integral thereby obtained formally is the first or- 
der density expansion of K[/] .) This is, indeed, the main 
issue of Sec. IIIID1 (For higher temperatures the QED 
scattering processes become dominant, and the complete 
treatment must be built on a quantum theory.) 

The infrared divergence reflects the long-ranged na- 
ture of electromagnetic interactions. The present case 
crucially differs from the nonrelativistic case in that the 
divergence exists also for the transverse electromagnetic 
interaction [see Eq. (|48p below] which is, however, only 
dynamically screened [7Ll4l|. As a result, the infrared di- 
vergence has to be healed by the short-ranged correlation, 
which is to be detailed in Sec. IIV Bl and was first noticed 
by Klimontovich in the formulation of a nonmanifestly 
covariant theory [38[ . The infrared cutoff is therefore set 



by inverse correlation radius, i.e., 

&min ~ Aq 1 . (45) 

The infrared divergence of the Landau collision integral 
roots in that the weak coupling approximation fails to 
capture the collective dielectric effects which leads to the 
short-ranged two-body correlation. This may be most 
readily seen by writing the relativistic Landau equation 
in the observer's frame. For this purpose, let us introduce 
the distribution function fo(x,p) defined as 

f(x,p) = f o (x,p)2m0( P o )6( P 2 - m 2 ) . (46) 

To simplify discussions we consider the spatially homo- 
geneous case. [The resulting / is denoted as fo(p,t) .] 
Then, the kinetic equation (fl7|) with the collision integral 
Ki [/] can be exactly rewritten as [38| 

| /o(p , t) = A.y„ p ' Q .(|._| 7 ) /o(p , i)A(p ', f) . 

(47) 

Here the tensor (in the 3-dimensional Euclidean space) 
Q reads out as (v and v' below are the velocity.) 

Q = 2e 4 J dujj d 3 k(5(cj-k-v)5(cj-k-v') 
kk f I [(k x v) ■ (k x v')] 2 \ 

X |k| 4 t| £ H(^,k)| 2 + | W 2 £^,k)-k 2 | 2 J ' ( ' 

where e"(w,k) = e ± (uj,'k) = 1 . The second term in the 
curl bracket arises from the interaction mediated by the 
transverse electromagnetic field. It is negligibly small 
in the nonrelativistic limit, i.e., |v| , |v'| <C f , while the 
first term arises from the interaction mediated by the 
longitudinal electromagnetic field and possesses the gen- 
eral structure of the Landau-Balescu-Lenard equation 
[20L [39l | . In comparing with the nonmanifestly kinetic 
equation obtained by Klimontivich [38| Eq. (|4"5)) suggests 
that the longitudinal (transverse) permittivity e" (e^) is 
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1. Formal collision integral 
We first introduce the following propagator: 

oo 

U ijz = \ n (ij\CM°(z)[2'Clh (z)} n \ij) . (50) 



FIG. 2: Diagrams leading to the Boltzmann collision integral. 



unity and, thus, the collective dielectric response is com- 
pletely ignored. Nevertheless this is an artifact of the 
weak coupling approximation. In fact, the nonmanifcstly 
covariant theory [381 ] has shown 



e"(w ->0,k) - 1 oc 



e x (w -> 0,k) - 1 oc 



1 



k|A D ) 2 

2 



ie 



(49) 



which suggests that the fluctuations of the longitudinal 
(transverse) electromagnetic field acquire a static (dy- 
namical) mass oc e 2 (namely the interaction strength). 



By inserting Eq. (|4*9")) into Eq. (|4*5)l it is clear that for a 
collision integral to account for the collective dielectric 
effects to go beyond the weak coupling approximation is 
inevitable. More precisely, (in the Newtonian physics) 
the collective dielectric effect is well known to result in 
a short-ranged two-body correlation, and to be responsi- 
ble for by more complicated correlation associated with 
higher order terms in the density expansion of the colli- 
sion integral [2pj . 

In the present context to treat the infrared divergence 
accurately we need to sum up all the so-called ring di- 
agrams (to be defined in Sec. IIVB[) . In Fig. Q] (b) we 
present the lowest order ring diagrams. In doing so we 
expect to obtain a manifestly covariant generalization of 
the Balescu-Lenard equation, where the infrared diver- 
gence is healed by collective dielectric effects. However, 
to carry out this program is far beyond the scope of the 
present work. We here limit ourselves to confirm this 
important observation in the case of global equilibrium, 
which is detailed in Sec. II V Bl 



D. Two-body scattering approximation: 
Relativistic Boltzmann equation 



In this part we wish to justify the relativistic Boltz- 
mann equation. The derivation below may be generalized 
to low density systems with short-ranged interactions, 
where the two-body scattering dominates. 



Upon passing to the time representation: Uij Z — > Uij (s) , 
we find 

^-U l3 (s) = (4° + £° + A4-) Uij (s) , Vki (0) = 1 ■ (51) 

Then, for Eq. (|33f we keep the entire A-expansion (Fig. [2]) 
which is called the two-body scattering approximation. 
The collision integral thereby obtained is denoted as 
IK2 [/] which reads out as 



X ds dS^d p 2 u 2 

Jo Jv 2 ®u% 

x£' 12 Ui2(s)C' 12 \Jl 2 (-s)f(x 2) p 2 )f(xi,pi 



(52) 



It differs from Eq. (|40|) in that the propagator U5 2 (s) 
sandwiched by £' 12 is renormalized into Ui 2 (s) (Fig. [21 
lower panel). 

The derivation is exact so far. In the leading order 
hydrodynamic expansion the spacetime inhomogeneity of 
f(x,p) and two-body scattering are decoupled. That is, 
we may also apply Eq. (pJ9")> to the collision integral K 2 [/] 
and, furthermore, approximate the operator U5 2 (— s) in 
K 2 [/] by unity. As a result, with s integrated out Eq. 
is simplified as 



»[/] = A 2 / 



dEfiid p 2 u 2 



(53) 



x£' 12 ZC' 12 f(x 1 ,p 2 )f(x 1 ,p 1 ) , 

where Z = J °° dslA\ 2 {s) . In order to find Z we use 
Eq. ([5Tj) to set up the following equation: 



Z = G + G \C' 12 Z , 



G 



dsU° 2 (s) 



(54) 



2. Two-body scattering and collision integral 

As before we employ the relativistic impulse approx- 
imation. There we insert Eq. (|34p into the two-body 
interacting Liouvillian A£' 12 . To proceed further we no- 
tice that to describe the free motion of particles 1 and 
2 , with the 4-momcntum vector p\ and p 2 , respectively, 
we may work in the center-of-momentum (CM) frame 
(Correspondingly we use the prime to denote vectors in 
this frame.) where p' : + p' 2 = 0. Because the corre- 
lation radius is short-ranged they do not interact with 
each other until at a distance of ~ Ad , when they start 
to repel each other. Eventually as their deviation reaches 
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order of Ad they undergo free flight again. Then, we as- 
sume that in the CM frame the interaction between two 
particles is switched on simultaneously. This, indeed, is 
perfectly legitimate because the hydrodynamic approxi- 
mation washes out effects arising from the small (coor- 
dinate) time mismatch ~ Ad of two particles. (Such an 
assumption may be released if we simplify short-ranged 
interactions as point-like collisions. In the latter case it 
can be shown that the derivations below are exact.) As 
such we find that, in the CM frame, 

Pi(O+P2(*') = (55) 

for all t' , moreover, the two-body interacting Liouvillian 
does not depend on the coordinate time t' any more: 

\C' t] {x[,p[- x'^) = A£^(xJ - x;,p^) . (56) 

That is, in such frame the two-body dynamics may be 
reduced into the dynamics of single particle subject to 
some external time-independent potential. The detailed 
analysis is presented in Appendix I A 31 

Let us introduce the following notations: v! = Pi/m = 
(u'°,u'), r> = (t',x',y' : z') = (t[ = t' 2 ,( x [ - x' 2 )/2) as- 
suming that the particle 1 approaches to (departures 
from) particle 2 from z' — —oo (at z' — +oo). Since XC 12 
is time-independent in the CM frame, from Eq. (|53|) we 
see that in the leading order hydrodynamic expansion the 
derivative with respect to it involved in the definition of 
G , namely Eq. (|54|) may be dropped out. As a result, the 
matrix element of G , denoted as G(r' ,p'; f' ,p') satisfies 

2u' • ^ G(r',p';r',p') = 6&(r' - v')5^\p'-p') . (57) 

Using this equation it can be checked that the matrix 
element of Z , denoted as Z(r',p'; r',p r ) solves 

|2u'- ^-X£' 1 ^Z(v',p';r',p') 

= S^(r' -v')8^\p' -ft). (58) 

Analogous to quantum mechanics Eq. (|58p describe the 
scattering of an incident "wave function" if- ln {p') under 
the "potential" which now reads out as — A£' 12 . And the 
out-going "wave function" tp out (r' ,p') is 

VWt(rV) =^ m (p') (59) 

+A J d 3 ? I dVu'°Z(r',p';f',p')A 2 ^(p')- 

I 



Using Eq. (J58]) we find 

|2u'~_A4 2 |v> OT t(r' 1 j/) = 1 (60) 

In deriving this equation we keep in mind that the two- 
body scattering is a local event with a characteristic 
scale Ad ■ Over such a scale incident wave functions are 
strongly scattered (by the potential — A£' 12 ) and, as such, 
the out-going wave function generally acquires a strong 
dependence on r' . In contrast, as the two-body scat- 
tering (in the CM frame) concerned the incident wave 
function -0 in in Eq. (|59|) may be regarded as a spacetime 
independent object. Nevertheless we notice that this is 
not true in the observer's frame where, instead, the inci- 
dent wave function reads out as 

^inO') = f{xi,Pi)f(xi,p 2 ) (61) 

varying over a macroscopic scale > Ad • In Eq. ([BTjl p' is 
uniquely determined by p\$ . 

Substituting Eq. (|59| into Eq. ([53]) we find that, in 
terms of the CM frame coordinates, the collision integral 
is 

Mf] = *j <?V' j u4 d 4 p'V° l' 12 {i>ont{v' \ P ') ~ ^in(p')} 

= xj dV^d 4 pV°4 2 Vw(rV) 

= J d*j/ Jd s r'2u'°n'~^ t (r',p'), (62) 

where in deriving the second line we use the fact that 
/ d 3 r'£' 12 — > an d in deriving the third line we use 
Eq. 



3. Boltzmann collision integral 

For Eq. (|62"j) let us integrate out r' first. For this pur- 
pose we fix the z'-axis to be in the direction of u' . Then, 



J d 3 r'2 M 'V • A^rV) = J dV2 W '°|u'| ^out(rV) 

= J dx'dy'{u'°\u'\ - u'°(-\u'\)} {^ out (r',p%^+oo ~ ^o»t(r',P% 



(63) 
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Remarkably, Eq. (|63|) depends only on the wave functions 
at z' — ±00 and is irrespective of the details of the two- 
body scattering. Noticing that the interaction vanishes 
at z' — > ±00 we thus find V , out(r',p')lz'->-°° = i>in(p') ■ 
In the CM frame applying the Liouville theorem further 
gives 

^out(r'y)U'-+oo = V^nW)) 

= f{x l ,p l )f{x u p 2 ). (64) 

Here p'{Sl') is the out-going 4- momentum vector depend- 
ing on the scattering angle f2' . In deriving the second 
equality we use Eq. (fBTj) and pass to the observer's frame 
with p\ t 2 standing for the out-going 4- momentum vector 
of the two particles. To proceed further we introduce the 
differential cross section in the CM frame er cm defined as 
dx'dy' = a cm (p\,p2 — > pi, p2)dQ,' . By definition, <7 cm is a 
manifestly covariant concept. Moreover, in the observer's 
frame the quantity u'°|u'| — u'°(— |u'|) may be written as 
\/ {ui ■ U2) 2 — 1 . Collecting everything together and tak- 
ing into account Eq. (|46|) . we obtain 

I - / rfS M2 rf 4 P2<A£' 12 / (2) 1 / (1) 

= K a [/o]. (65) 
Here the collision integral is 

Mf} = J J dn'^a cnl ^/( Ul -U2) 2 -l (66) 

x{fo(xi,pi)fo(xi,p 2 ) - f Q (x 1 ,p 1 )f (x 1 ,p 2 )}. 

Eq. (|66p is the well known relativistic Boltzmann colli- 
sion integral [l(| EU , where the details of the two-body 



scattering enter through the invariant differential cross 
section u cm . 



IV. CORRELATION AT RELATIVISTIC 
MANY-BODY EQUILIBRIUM 

It is important that the perceptions of a macroscopic 
observer correspond to the hierarchy of physical correla- 
tion functions namely 91 oo rather than D , and the hier- 
archy is fully determined by the solution to the general 
transport equation (jTTJ) with the collision integral given 
by Eq. (fl"8|) . Obviously, the processes with regard to the 
physical distribution described by Eq. (I17|) are observer- 
independent which, in the chosen inertial frame, may be 
parametrized by the coordinate time. For rarefied plas- 
mas Eq. (fT"7|) becomes a local equation, and the collision 
integral generally admits the so-called "collision invari- 
ants" which, typically, are <fi(p) = l,p M , i.e., 

/ d i p^(p)K[f(x,p)}=0. (67) 

They amount to the particle number and the energy- 
momentum conservation law, from which macroscopic 
relativistic hydrodynamics follows. 

As the system reaches global or local equilibrium the 
physical one-body distribution function nullifies the colli- 
sion integral, i.e., K[/ cq ] = . The former differs from the 
latter in the spacetime independence of / cq . By inserting 
/oq into Eq. (fl~9)) we obtain 



(r?(V ' ' >i)|9*°o = / ds H / dE^Piuf T\f e q(x k ,Pk) 

x (1^(1, ■ • • , i)|C£t°(s - S ')£0OVexp(-sVTV)|l| • • • |n) 
I 



(68) 



which, in principle, gives the entire hierarchy of physical 
equilibrium correlation functions. Therefore, the conjec- 
ture of Hakim on the many-body equilibrium @ is jus- 
tified. In Ref. 9] only the first component, i.e., (12|9Too 
is formulated using the weak coupling approximation. 
This hierarchy, together with f eq (x,p) fully captures the 
many-body equilibrium of rarified relativistic plasmas. It 
allows one to go beyond the kinetic approximation 0, Q 
to formulate the relativistic many-body equilibrium bas- 
ing on the action-at-a-distance formalism. 



A. Jiittner equilibrium 



Let us now exemplify these general principles in the 
case of 1Ki ; 2 [/] , which can be easily shown to pre- 
serve Eq. (|67[) . This then results in the relativistic 
hydrodynamics-a well established subject [Til ] which we 
shall not proceed further. Encompassed by the hydro- 
dynamics the system irreversibly evolves into the local 
Jiittner equilibrium. Indeed, define the local entropy flux 
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S fl (x) as follows: 



S^x) = - / d 4 p2me{p )5{p 2 -m 2 )u' i f (x,p)hxf (x,p), 

(69) 

where fo(x,p) is defined by Eq. (|46[) . Then, from the 
collision integral the iJ-theorem follows JJJ, i.e., 



tr a (a!) = fl M 5"(as) >0. 



(70) 



The local entropy production a s (x) vanishes wherever 
fo{x,p) reaches the local Jiittner equilibrium: 



fj(x,p) 



p{x)(i{x) 



Airm 2 K 2 (m[3(x)) 



(71) 



Here p(x) is the invariant particle number density, K 2 is 
the modified Bessel function of order two, and /3 M (x) is 
the (local) reciprocal temperature timelike 4- vector defin- 
ing the (local) temperature T{x) through l/k B T(x) = 
(3(x) 2 . Suppressing all the spacetime dependence of 
the parameters of fj(x,p) further leads to the (global) 
Jiittner equilibrium: 



/e q (p) = fA P )2m9(p°)S(p 2 - m 2 ) 



(72) 



The debate on the special relativistic version of 
the Maxwell-Boltzmann distribution has stemmed from 
Ref. [l7t | where a different Boltzmann equation, though 
manifestly covariant, is proposed. There, proceeding 
along the line similar to Eqs. ([69 |) -([7T j) results in an alter- 
native one-body equilibrium distribution. In the present 
work by the proof of the relativistic Landau and Boltz- 
mann equation namely Eq. (| 1T[) with the collision in- 
tegral ([4"Tj) and (|66|) is implied the following important 
fact: In rarefied plasmas the Jiittner distribution is no 
longer a phenomenological hypothesis, rather, is well jus- 
tified at the full microscopic level and, importantly, suits 
the covariant principle. Thus, we provide a solid sup- 
port for very recent 1-dimensional numerical simulation 
[26j and advocates the critical analysis on alternatives 
to the Jiittner distribution [13] • Notice that the spatial 
dimensionality in the numerical simulation differs which, 
we believe, is of minor importance. (In fact, for plas- 
mas the relativistic many-body dynamics and thereby 



the collision integral are substantially simplified in the 
1-dimensional geometry. There, the electromagnetic in- 
teraction is purely longitudinal.) What is crucial is that, 
there, the underlying many-body dynamics is dictated by 
point-like collisions and is well described by the action- 
at-a-distance formalism and, as a result, the entire scope 
of the present manifestly covariant correlation dynamics 
applies. 



B. Equal-time two-body correlation 

We further proceed to calculate the correlation func- 
tion at the Jiittner equilibrium. For simplicity let us con- 
sider the simplest two-body correlation function, which is 
C eq (xi,pi;x 2l p 2 ) = (12|91oo . Here we use the subscript 
"eq" to denote that the physical one-body distribution 
is at global Jiittner equilibrium, i.e., f(xi,pi) — / oq (pi) • 
In particular, to investigate the behavior of the spacelike 
correlation function we will calculate below the so-called 
equal-time correlation function. 



1. Bare correlation function 

To the lowest order interaction expansion [see Fig. [3] 
(a)] Eq. (|68| gives a bare two-body correlation function: 



CeqOl. Pi \X2,P2) 



ds (12|U°(s) A£'|l|2) / eq (pi)/ eq (P2) • (73) 



Under the relativistic impulse approximation (see Ap- 
pendix |X2] for details) it is translationally invariant, i.e., 
C cq (x 1 ,p 1 ;x 2 ,p 2 ) = C° q (a;i - x 2 ,Pi,p 2 ) ■ Indeed, insert- 
ing the matrix elements of ilo an d A£' into Eq. (173[) we 
find 



d*k 



C£,(n - x 2 , PuP2 ) = I e *(*i-* 2 ) c° cq (k, Pl , P2 ) 



(74) 



with 



z87r 2 e 2 1 

C°(k,p!,p 2 ) = T2~~i — 1 ?{5(k ■ Pz)[k ■ (3(pi ■ P2) - (3 ■ p 2 (k ■ pi)] 

H k z ik ■ (pi — p 2 ) 

-S(k ■ pi)[k ■ f3( Pl ■ pa) - p ■ Pl (k ■ P2)]}/eq(Pl)/eq(p3) , (75) 

where the wave vector fc M = (uj, k) . 

Since the equilibrium is global we may take the advantage of the Lorentz invariance of Eq. (|74p and choose the 
observer's frame where /3 M = (l/fc^T, 0) , and the spacelike 3-surface is the usual 3-dimensional Euclidean space. In 
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this frame the bare equal-time correlation function is 

C c q ( Xl -x 2 , Pl ,p 2 ) = J * e^to-*,)^ 

C° q (k, Pl ,p 2 ) = J ^C c >,k, Pl ,p 2 ). 
Notice that in the observer's frame Eq. ([75]) is written as 



Pl,P2) , 



1&7T C 



5(k ■ p 2 ) 



+S(k-pi) 



u (Pi ■ P2) 

fc B Ti(wp?-k-pi) ^ k B T 

u) (pi • p 2 ) ipi 



Integrating out to we find 



Co q (k,pi,_p 2 ) oc 



k B Ti{cjp° 2 -kp 2 ) k B T 

/cqbl)/eq(P2) 



(76) 



/cqbl)/cq(P2) • 



(77) 



(78) 



where the irrelevant numerical overall factor is ignored. Eq. (|78p shows that the bare equal-time two-body correlation 
is long-ranged, i.e., C° q (xi — x 2 ,pi,p 2 ) oc |xi — x 2 | _1 reflecting the long-ranged nature of electromagnetic interactions. 

2. Short-ranged equal-time correlation: Collective dielectric effects 



To go beyond the weak coupling approximation we extend the so-called ring approximation [34| to the relativistic 
context. That is, we will sum up all the so-called ring diagrams. [Typical diagrams are given in Fig. [3] (b)-(d).] 
They are obtained from Fig. [3] (a) in the following way: The propagating lines to the left of the vertex-associated 
with no particle annihilation-are dressed by a sequence of vertexes. For each of them an additional particle joins the 
vertex from the right and continues propagating to the left, and eventually is annihilated at the next vertex. (It can 
be shown that the propagating lines to the right suffer no renormalization effects.) In selecting these diagrams it is 
implied that in Eq. (|68p exp(— sVTV) is set to be unity. Summing up all these diagrams we find that the two-body 
correlation function solves the following Dyson equation: 



C eq (x 1 ,p 1 ;x 2 ,p 2 ) = C° q (x 1 , Pl ;x 2 ,p 2 )+ / ds e' s{ < d ^ + < 9 ^ 



£ 3 <gi[/f 

x {^'i3fcq{Pi)C cq {x 3 ,p 3 - 1 x 2 ,p 2 ) + A£ 23 / cq (p 2 )C oq (x3,p 3 ;a;i,pi)| . (79) 

Let us now pass to the observer's frame chosen above, where t\ = t 2 = t 3 and calculate the equal-time corre- 
lation using Eq. J79|) . At global equilibrium such a correlation function does not depend on the coordinate time. 
Consequently, upon passing to the spatial Fourier transformation we obtain 

S(k- P3 ) 



/duj 1 f 
ji z jkp J d 



4 



ik • (pi 



P2J 



UJ 

ke^T 

LU 



(Pi • Pa) 



(P2 ■ pa) 



k B T 

_pL 

k R T 



(k-pi) 



(k- P2 ) 



/ eq (pi)C eq (k,p 3 ,_p 2 ) 



feq(P2)C e q(KP3,Pl, 



(80) 



where C cq (k,pi,pj) is the Fourier transformation of C cq (xi — x.j,pi,pj) . To further proceed we choose the direction of 
k as the a;-axis and integrate out u> . It is natural to expect that C eq (\i,pi,pj) possesses the spherical symmetry with 
respect to Pj j . As a result, we obtain 



C eq (k,pi,p 2 ) = C° q {k,p 1 ,p 2 ) 



47r 2 e 2 1 



1 



k B T |k| 2 pi x -p 2x 



Plxfeqipi) / d 4 p 3 M 3 C eq (k,p 3 ,p 2 ) 



P2xfcq(P2) / d 4 p 3 ulC eq (k,p 3i p 1 ) 



(81) 



r 



Inserting Eq. ([75)1 into it we find the solution to be 

/eq(Pl)/eq(P2) 



Coq(k,Pl,_p 2 ) OC 



|k| 2 + A D ^ 



(82) 



Eq. (|82|) with Ad given by Eq. (|28|) fully agrees with 
the result of Klimontovich who used a completely differ- 
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FIG. 3: Diagram of bare two-body correlation (a) and typical 
diagrams leading to short-ranged correlation (b)-(d). Solid 
circles stand for the Jiittner equilibrium. 

ent (nonmanifestly covariant) approach [381 ]. They sug- 
gest that the equaltime two-body correlation function (in 
the phase space) displays an exponential decay with the 
correlation radius Ad , which is none but the nonrela- 
tivistic Debye length. The relativistic insensibility (rest 
mass) of Ad was noticed long time ago pfjl |. Comparing 
Eq. (|82p with Eq. (|78|) we find that despite of the appear- 
ance of the transverse electromagnetic interaction, which 
is screened in a dynamical manner 0, EH ] , the collective 
screening of the longitudinal electromagnetic interaction 
renders the equal-time correlation function short-ranged. 
Moreover, the static longitudinal permittivity is found to 
be 

^-°>v- 1 = wkr- (83) 

Thus, the ring diagrams suffice to describe the collec- 
tive dielectric effects, justifying the introduced lower cut- 
off Aq 1 in the relativistic Landau collision integral Ki [/] . 
Of course, to heal the infrared divergence of Ki [/] accu- 
rately we need to sum up all the nonequilibrium ring 
diagrams as depicted in Fig.Q](b). This is technically far 
beyond the scope of the present paper and we leave for 
future studies. 



V. CONCLUSIONS 

The present work is motivated by very recent physical 
and mathematical progresses in 1-dimensional relativis- 
tic many-body systems 0, |42j . Although the attempts 
of reconciling the correlation dynamics and the special 
relativity were undertaken by the Brussel- Austin school 
long time ago, the relativistic correlation dynamics was 
formulated in a nonmanifestly covariant manner [12l . 1 1 61 ] . 
Because of this drawback such a theory has been less ad- 
vanced so far with regard to the practical applicability. 
The possibility of formulating a manifestly covariant cor- 
relation dynamics was first discussed by Israel and Kan- 
drup [2|, |3j ■ Unfortunately, there have been no progresses 



toward this direction. In this series of papers we sub- 
stantially extend the analysis of the relativistic classical 
nonequilibrium statistical mechanics of full interacting 
many-body systems by Israel and Kandrup in the con- 
text of special relativity. Then, a manifestly covariant 
correlation dynamics results, which bridges macroscopic 
phenomena such as hydrodynamics and equilibrium and 
microscopic deterministic relativistic many-body dynam- 
ics. It is a statistical theory of ensembles of a bundle of 
world lines, and is suitable for studies of dynamics of dis- 
tribution functions of full classical many-body systems. 
For a summary, we carry out the following program: 

(i) For a classical interacting system composed of N 
particles the well established action-at-a-distance for- 
malism naturally introduces an 8-ZV-dimensional T phase 
space, for which we introduce a probability distribution 
function. From the conservation law a single-time Liou- 
ville equation follows, describing the manifestly covari- 
ant global evolution of the distribution function of the T 
phase space. 

(ii) For the manifestly covariant single-time Liouville 
equation we perform the correlation dynamics analysis. 
First, we introduce the correlation pattern representa- 
tion, where the cluster expansion of reduced distribution 
functions is formulated in a manifestly covariant manner. 
Then, it is found that the evolution of the full A-body 
distribution function (of the T phase space) may be re- 
duced into that of the one-body distribution function (of 
the 8-dimensional fj, phase space) which is described by 
a closed nonlinear equation. However, the evolution of 
such a one-body distribution function may not be ob- 
servable. Rather, the stationary solution of this closed 
nonlinear equation results in an exact general transport 
equation. The solution of the latter is physically observ- 
able and fully responsible for macroscopic hydrodynam- 
ics and, very importantly, determines the entire physical 
correlation functions. Such picture is encapsulated in 
Eqs. (fl~7f and (|19p . which constitute a manifestly covari- 
ant version of the Bogoliubov functional assumption-the 
general principle underlying kinetic processes of many- 
body systems. 

(iii) The collision integral of the general transport 
equation is expressed in terms of the density expansion. 
Therefore, for low-density systems only the leading term 
needs to be taken into account. For such obtained colli- 
sion integral the perturbative theory with respect to the 
interaction is formulated. Then, under the relativistic 
impulse approximation we recover the Landau collision 
integral by keeping the leading order interaction expan- 
sion, while recover the Boltzmann collision integral by 
keeping the entire interaction expansion. 

(iv) The analysis of the entropy production due to the 
Landau and the Boltzmann collision integral shows that 
a low-density system (for example, rarificd plasmas) tend 
to be driven to the (local) Jiittner equilibrium. Replacing 
the physical one-body distribution function in the hierar- 
chy (|19[) with the (global) Jiittner equilibrium, we justify 
the long standing many-body equilibrium conjecture of 
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Hakim. In particular, we calculate explicitly the (static) 
two-body correlation function for plasmas and analyze 
the collective screening effects. 

We remark that the emergence of a manifestly covari- 
ant formulation of the Bogoliubov functional assumption 
is by no means obvious. Indeed, in relativity the break- 
down of the simultaneity as a covariant concept pro- 
foundly affects the mathematical structure of the cor- 
relation dynamics. Among the most important modifi- 
cations is that the microscopic states, in general, are no 
longer the representation points in the T phase space, 
rather, are a bundle of N (segments of) particle world 
lines. (There do exist exceptions, for example, the case 
studied in Ref. ,42].) Indeed, the introduction of the T 
phase space is merely suggested by the form of the motion 
equations. In particular, that the force does not depend 
on the acceleration of the acted particle (Thus, the radi- 
ation reaction is ignored as electromagnetic interactions 
concerned.) is crucial to the present choice of the T phase 
space. As a result, the building block operator namely 
the two-body interacting Liouvillian and the distribution 
function are determined not only by the positions and the 
momenta of two participating particles, but also by their 
world lines (more precisely, the phase trajectories) pass- 
ing through the given positions and momenta. The very 
origin of the Bogoliubov functional assumption namely 
Eqs. (TTT)) and ([TO)) is that the deterministic many-body 
dynamics lose the memory of the given (segments of) 
particle world lines at large global proper times. 

In fact, it is the nonmanifestly covariant version of the 
Bogoliubov functional assumption, formulated in a rather 
straightforward manner, that is commonly adopted in 
the literature of plasmas [HI, H] ■ However, the proof 
has been lacking. The manifestly covariant Bogoliubov 
functional assumption namely Eqs. (TTTj) and (fTi?)) justi- 
fies such an formulation. Indeed, it was first noticed by 
Dirac, Fock, and Podolsky [45| and stressed by Israel Q 
that it is legitimate to identify the coordinate times of 
particles, i.e., to set t\ = tq. = ■ ■ ■ = t^ = t for the 
manifestly covariant equation £J\f = . In doing so we 
obtain Eqs. (JT7J) and (TT5)) with all the coordinate times 
set to be the same. Next, as (implicitly) done in various 
nonmanifestly covariant theories we again apply the rel- 
ativistic impulse approximation. The key point is that 
such an approximation does not destroy the (manifest) 
covariance of Lorentz force. Then, the collision integral 
(I33p remains unaffected except that there t\ and t2 have 
to be identified. For this collision integral let us keep only 
the leading order interaction expansion. Immediately, we 
arrive at a relativistic Landau collision integral which is 
a nonmanifestly covariant version of Ki[/] . Thus, we ex- 
plain the insightful observation made in Ref. 0], where 
it was pointed out that to derive relativistic transport 
equations by using a nonmanifestly covariant approach 
may not be a problem as electromagnetic interactions 
concerned. 

We stress that upon passing from the motion equations 
and @ to the single-time Louville equation © some 



information of the underlying many-body dynamics is 
lost. Thus, they are not equivalent. The manifestly 
covariant correlation dynamics stems from the latter. 
Surprisingly, for rarefied plasmas the present correlation 
dynamics recovers the Landau and the Boltzmann 
equation which are fully responsible for various macro- 
scopic phenomena such as hydrodynamics and the (local) 
equilibrium. Thus, we are led to the following conjecture: 

For systems with arbitrary macroscopic particle number 
density the manifestly covariant single-time Louville 
equation (0) completely determines macroscopic rela- 
tivistic hydrodynamics. 

Although the relativistic Landau equation has been 
justified before especially by the manifestly covari- 
ant nonequilibrium statistical mechanics developed in 
Refs. [U, y], compared to the earlier theories the mani- 
festly covariant correlation dynamics presented here has 
the additional advantage in going beyond the weak cou- 
pling approximation: First, it clearly shows that the rela- 
tivistic Landau equation is the leading order interaction 
expansion of the relativistic Boltzmann equation. It is 
the two-body scattering approximation that is used in 
deriving the latter equation, and the collision integral 
(with the accuracy of the leading order density expan- 
sion) is expressed in terms of the interaction expansion. 
The present theory allows one to sum up this expansion 
and, as a result, the scattering cross section is beyond the 
weak coupling level. Second, the general collision integral 
(TTS)) and the hierarchy of correlation functions (flU)) fur- 
ther allow one to systematically go beyond the two-body 
scattering approximation. This is exemplified in the cal- 
culations of the static two-body equilibrium correlation. 
It is promising in understanding the infrared divergence 
of the relativistic Landau collision integral where collec- 
tive dynamic screening phenomena play crucial roles, and 
in generalizing the Balescu-Lenard equation [2(| so as to 
suit relativity principles. (To our best knowledge, so far 
only the nonmanifestly covariant Balescu-Lenard equa- 
tion has been obtained [H, [H, 5(|.) Finally, the present 
theory allows one to explore the anomalous transport 
phenomenon arising from the long-time tail correlation 
effects [47] . These issues seem to be far beyond the reach 
of earlier theories such as Refs. 0, In- 
closing this paper we remark that although this series 
of papers focuses on the classical correlation dynamics, 
it may have important implications on the relativistic 
quantum nonequilibrium statistical mechanics. Indeed, 
to formulate the latter at the level of full many-body 
quantum dynamics has been of long term interests and 
remained controversial so far [48j] . Besides of the man- 
ifest/nonmanifest covariance issue another severe diffi- 
culty arises as the quantization of relativistic many-body 
classical dynamics concerned. In order not to lose the 
manifest covariance Schieve recently started from the so- 
called covariant Hamiltonian dynamics and formulated a 
theory of relativistic quantum statistical mechanics [48| . 
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It is important to notice that the quantum Boltzmann 
equation derived there is with regard to "events" in the 
spacetime rather than particles. For this reason nei- 
ther can this equation be justified by the microscopic 
approach [14| basing on the generalization of traditional 
nonequilibrium Green's function theories (49l . l5(il ] or by 
refined relativistic nonequilibrium quantum field theories 
[SHE], nor the classical limit of this equation be justified 
by other classical microscopic theories [1, H] ■ In contrast, 
the quantum transport equations obtained by complete 
quantum treatment, in the classical limit, agree with the 
transport equations shown in the present paper. This 
justifies the legitimacy of basing classical relativistic sta- 
tistical mechanics on the action-at-a-distance formalism. 
Many natural problems thereby arise: May it be possible 
to be extended to the quantum case? What is the quan- 
tum counterpart of the proper time-parametrized evolu- 
tion? In particular, over a field theory the advantage of 
proceeding along this line to formulate a quantum sta- 
tistical theory is to allow for nonperturbative treatment 
and to encompass some global properties, as admitted 
by many authors (for example, see Ref. [la]). Unfortu- 
nately, to solve these problems is by no means an easy 
task because of the lack of Hamiltonian of the action-at-a- 
distance formalism. A possible prescription is to enlarge 
the configuration space to accommodate the field degrees 
of freedom. Treating the particles and the fields on the 
same footing one may proceed to construct a many-body 
Green function by adopting the prescription of Feynman 
[EH H2| j namely to express it in terms of the functional 
integral over all the paths of both particles and fields 
parametrized by separate proper times. Furthermore, it 
is expected that such a many-body Green function satis- 
fies a Stiickelberg-type equation [53] suitable for the use 
of the Wigner function technique. As such the entire 
theoretical scope of nonrelativistic quantum correlation 
dynamics 34] might become applicable. We leave the 
detailed analysis for future studies. 
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APPENDIX A: RELATIVISTIC IMPULSE 
APPROXIMATION FOR ELECTROMAGNETIC 
INTERACTIONS 

In this appendix we wish to study the relativistic im- 
pulse approximation to two-body interacting Liouvillian 
for classical electromagnetic interactions. 



1. Exact two-body interacting Liouvillian 

We notice that, in general, the interacting force F£ has 
the form as follows [23, |54| : 

/+oo 
dT 3 s(p ij )J 7f "'\ Xj=Xj ( Tj ) p ui , (Al) 
-oo 

Pij — ( x fii ~ x nj)( x i ~ x j) j 

where is an antisymmetric tensor, and the role of 
function g(pij) is to invariantly connect Xi with one (or 
several) points at the world line Xj(jj) . Let us exemplify 
Eq. (|A1|) and derive the exact form of two-body inter- 
acting Liouvillian in classical electrodynamics. For this 
purpose we employ the Wheeler-Feynman formalism [32| 
and start from the Fokker action which crucially allows 
us to eliminate the self-action: 



S = -m^T / ds^u^u 1 *) 1 / 2 - e 2 ^ // ds t ds 3 u^u^ 5(p i3 ) (A2) 
i J i<j JJ 



with u^(s) = dx^{s)/ds. 

Then, we apply the variation principle to this action. Demanding 



SS = to V / ds~ $— 5x^ + e 2 V // dsidsj 5( Pij )] - 9f [tf 6( Pij )}} u^Sx^ = (A3) 

J d Sl ^u vi u\ f^JJ 



and defining dri = ^Ju^u^ dsi and = uf / u v iu" , we where F^ v is the electromagnetic tensor adjunct to par- 
obtain tide j , and all the field point is set to be the position 

Fg{xi) - [d»AXx)-d"A»(x)]\ x=Xi) (A4) 
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of particle i at , i.e., x = Xifc) . The vector potential, 
adjunct to particle j , is given by 

A$(x) = e y $((a; - x,) 2 ) u£ dr, . (A5) 

Notice that here both cc^ and depend on the proper 
time Tj . Eqs. (|A4[) and (|A5[) describe the classical elec- 
trodynamics of electronic systems in terms of direct in- 
terparticle action. Notice that the vector potential given 
by Eq. (|A5|) satisfies Lorentz gauge: 

d M A!j(x) = 0. (A6) 



Inserting Eq. into Eq. (|A"4|) gives Eq. ([A"Tj) with 

J*» = - xfW - tftf - x))] , (A7) 

s(Pij) = S'(Pij) • (A8) 

Here S'(x) = dS(x)/dx . 

With the help of Eqs. ([A"4"|) and (fX5j) the exact two- 
body interacting Liouvillian is written as 



XC' ij (xi,p i ;x j ,p j ) = -e j dr in^uf (r)d? - u^(T)d^}S(( Xi - x {t)) 2 ) 

+u vj [uUr)df - <(r)9J] 5(( Xj - x t {r) f) A 



(A9) 



2. Two-body interacting Liouvillian: relativistic impulse approximation 

Now we proceed to find the relativistic impulse approximation to the two-body interacting Liouvillian. For this 
purpose we use the Dirac identity: d fJl d fl S(x tl x fl ) = 47r<5( 4 )(a;) to obtain 



d A k e tkK 
(2^)4 "fe2 



8{ x ^) = -^l— T _. (A10) 



Substituting it into Eq. (|A9j) and integrating out r , with Eq. (|34|) taken into account we obtain 

F$(x i ,p i ) = — j — ^ p *(*-Pi)[^(Pi-Pi)-J#(fc-Pi)] (AH) 

I 



and likewise for F^(xj,pj) . Eventually we find the two- 
body interacting Liouvillian, under the relativistic im- 
pulse approximation to be Eq. ([55]) . (Notice that there 
we need to make the replacement: 1 — ^ ?' , 2 — * j' for the 
subscript.) It is important to notice that the mass-shell 
constraint is preserved by the forces: 

Fij ■ i>, = Fji -pj = 0. (A12) 



3. Two-body scattering in the CM frame 

Suppose that the coordinate times of particle i and 
j are identified, i.e., ti — tj = t. Then, in order for 
Eq. (f54|) to accommodate a solution in the CM frame we 
need to establish the following (7^ = ^/m 2 + \ 2 /m) 

Lemma. For the relativistic two-body dynamics de- 



scribed by the following equations: 

dxi p l dpi 
at m at 

dt _ to ' ld dt ~ 31 [ ' 

with the force given by Eq. \A1 if at t — exists 
Pi(0) + Pj(0) = 0, then it holds for all t > 0, i.e., 
Pi(t)+Pi(t) = 0. 

Proof. This is quite obvious. First of all, enforc- 
ing ti = tj = t we find, from Eq. (jAllj) . that the 
force arising from interactions depends on x, — Xj and 
Pij . Then, Eqs. (|A13|) and (|A14p uniquely determine 
the phase trajectories: Xij(i) provided that the initial 
condition: Xjj(0) , Pij(0) is set. 

On the other hand, consider the trajectories satisfying 
Xi(t) - x = -(xj(t) - x ) = x(i) , whe re x = (x t (0) + 
x i (0))/2 . We find that they solve Eqs. (fAT3| and (|A14|) 
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by noticing 

87r 2 e 2 [ d 3 k sin(2k-x) 

" ~^~J\ k \>x^ (2tt)4 (k.p) 2 -|k| 2 £ 2 

x£[k(£; 2 + |p| 2 )-2(k.p)p], (A15) 

where p(t) = p,(i) = — Pj(i) satisfying p = 71EX 
(7 = \Jm 2 + |p| 2 /m), and E = \/m 2 + |p| 2 . Thus, 
Eqs. (IA13[) and (IA14P do admit the solution satisfy- 
ing Pi(t) + Pj(t) — irrespective of *-i,j(Q) . With the 
uniqueness taken into account the lemma follows. Q.E.D. 

Remark. The lemma shows that under the relativistic 
impulse approximation the two-body scattering may be 
described in the usual CM frame where pi (t) + pj (t) = 



and Xj (t) + Xj (t) = . Importantly, in this frame the 
two-body dynamics may be reduced into the dynamics 
of single particle subject to the external force given by 
Eq. (TA151) . 



APPENDIX B: RELATIVISTIC LANDAU 
COLLISION INTEGRAL 

We now calculate Eq. (|40ll under the relativistic im- 
pulse approximation. For this purpose we insert the 
two-body interacting Liouvillian namely Eq. (|36[) into it. 
Notice that XC' i j(xi,pi;Xj,pj) depends on Xij through 
Xi — Xj . Then, 



V m J Jo Jz 2 <s>m 



lo J^xii " *JJ (27r)M27r) 4 

ln 2 e 2 \ 2 f ,„ „ „ ff d 4 k d 4 k' 



dT lfl2 d i p 2 pi 2 l J J ■ 



xe ifc.( ai -« 3 ) g u{k) e ik'-( Xx - X2 ) 5{k , . pi _ k >. p2 ) g 12{ k>) f(x 1 ,p 2 )f(x 1 ,p 1 ) , (Bl) 

where in the derivation of the second equality we have integrated out s. Inserting the expression of Gx 2 (k') [see 
Eq. (f36|) ] into it we obtain 

Kl[/1 = .J^f) [ tzMllgL^L 

\ m J J^ m A J J (2tt) 4 (27r) 4 

xe «.(x 1 -x a ) g 12{k)e ik ■(«.-»,) S (k> ■p 1 )S(k'-p 2 ) (p x -p 2 ) —I^L-^j f(x 1 , P2 )f(x uPl ) . (B2) 

Let us carry out the spatial integral first. For this purpose we enjoy the manifest covariance and choose T, 2 to be the 
usual 3-dimensional Euclidean space. As a result, we obtain 



dY lll2 e l{k+k ' Hxi - x ^ ) = (27r) 3 (5 (3) (k + k') e ^+^')^~^) . (B3) 

where we keep in mind that the wave vector = (u>, k) . Then, 

I dYj^ 2 p£ e WM*i-*2) 6 ( k . pi) S (tf . pi) fitf . p2) 

dE M2 pZ e WY(xi-x2) S ( k . p2) . pi) g(tf . p2 ) = (2?r) 3 ^(4) (fc + ^ . pi )S(k-p 2 ), (B4) 
dE M2 fc" e <(*+*')-(«i-^) S(k ■ Pl ) 6(k' ■ px) S(k' ■ p 2 ) = . (B5) 



Let us substitute the expression of Qx 2 (k) into Eq. (|B2|) and take Eqs. (|B4|) and (|B5[) into account. Noticing 

[*(* • P2)W{px ■ p 2 ), d/drt] = [S(k ■ px)k"(px ■ P2 ),d/dp 2 l ] = , (B6) 

we arrive at Eq. (|4"1"T) . 
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APPENDIX C: DIVERGENCE OF TENSOR e M " 



Observing Eq. (|42|) it is easy to see that possesses the following structure 



(CI) 



enforced by the Lorentz invariance, with 01,2,3 the coefficients which are invariant under the Lorentz transformation. 
Then, from the identity: pi^e^ = y 2 ^ v — we find C\ = • u 2 ) 2 — l]c 2 , c 3 = — (u\ ■ u 2 )c 2 . Substituting these 
two relations into Eq. (jClj) gives 

(C2) 
(C3) 



= c 2 {[{ Ul ■ u 2 f - l]g" u + «< + u%ul) - (u x • m 2 )(« + «)} ; 
where C2 is determined by 

1, 



C2 = ^[(ui-u 2 ) 2 -l]- 1 g^. 



We now come to calculate ci . For this purpose we insert Eq. 



C2 



(mi • u 2 f 



1 



d k5(k ■ u\)5(k ■ 112) 



into Eq. (|C3[) arriving at 
1 



jfe • k 



(C4) 



To proceed further we enjoy the Lorentz invariance and choose the direction of Ui as the rr-axis of the 3-dimensional 
Euclidean space. With u> integrated out we find 



f'2 



(ttl • U 2 f 



7172 [("1 ' U2) 2 - 1] 
4 {P1-P2) 2 



7i72[(«i • U2) 2 - 1] \vi - v 2x \ 



.3i &{kx{vi - V2x) - k y V2y - k z V 2z ) 
{k x V X f - k2 



d 2 k t 



k±-V21 



7l(Dl-W2x) 



(C5) 



where in the last line we use the notation: = (k v , k z ) , V2_l = {v2y, V2z) ■ 

The integral in the last line may be easily performed by passing to the polar coordinate system, i.e., , 9) with 
the axis chosen to be the direction of V2j_ ■ As such the integral is factorized into the integral over 9 which is finite 
[Hll], and the integral over k± which suffers from both the ultraviolet and the infrared divergence. More precisely, 



d 2 k t 



k±-v 2 ± 



i(vi- 



dk± 
k± Jo 



dd 



1 



2tt 



y'l(v 1 -v 2 o : ) 2 



7f(t)l-t)2x) 2 

dkj_ 



(C6) 



We see that, indeed, a logarithmic divergence results. Substituting it into Eq. (|C5|) gives 



c 2 = —2ire 



-2Tie 4 



dk^ (ui ■ u 2 ) 2 
k± 7172 • u 2 ) 2 
dk^ (ui ■ u 2 ) 2 



1 



1 



k± 7l72[("l ' u 2 ) 2 - 1] ^-(1 - v()(\ - V 2 ) + (1 - Vl V2 x f 

I 



(C7) 



Taking into account the identity: u\-u 2 — 7172(1— viV2x) 
eventually we find 



c 2 = — 2ire 



4 f dk ± (ui • u 2 ) 2 



k± [(ui • u 2 ) 2 - I] 3 / 2 



(C8) 



Inserting it into Eq. (|C2[) we arrive at Eq. (T53]). 
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